Efficient C-Phase gate for single-spin qubits in quantum dots 
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Two-qubit interactions are at the heart of quantum information processing. For single-spin qubits 
in semiconductor quantum dots, the exchange gate has always been considered the natural two-qubit 
gate. The recent integration of magnetic field or g-factor gradients in coupled quantum dot systems 
allows for a one-step, robust realization of the controlled phase (C-Phase) gate instead. We analyze 
the C-Phase gate durations and fidelities that can be obtained under realistic conditions, including 
the effects of charge and nuclear field fluctuations, and find gate error probabilities of below 1CP 4 , 
possibly allowing fault-tolerant quantum computation. 

PACS numbers: 03.67.Bg, 03.67.Lx, 73. 21. La 



The spin of a single electron is a natural degree of free- 
dom for storing quantum information. In semiconductor 
quantum dots, all operations required for quantum com- 
putation with single spins have been demonstrated ex- 
perimentally in recent years: single-shot read-out [l[ and 
coherent rotations of individual spins Q as well as co- 
herent interactions between two neighbouring spins Q. 
Single-spin relaxation times, Ti, are extremely long, up 
to 1 s p], 0] • Recent measurements in GaAs quantum 
dots gave spin-echo coherence times, T2.ec/10, of up to a 
few hundred fxs |5j . In quantum dots built from materials 
without nuclear spins, even longer spin coherence times 
are expected. 

Despite this impressive progress, there have been no 
experimental reports of elementary two-qubit gates such 
as the controlled-NOT gate or the C-Phase gate. These 
gates are important since they are, together with single- 
qubit rotations, the natural building blocks for quantum 
algorithms. Implementing the quantum Fourier trans- 
form for instance, comes down to concatenating a se- 
quence of two-qubit C-Phase gates @ . 

In their seminal paper, Loss and DiVincenzo present 
a scheme to achieve a two-qubit C-Phase gate by two pe- 
riods of coherent evolution under the exchange interac- 
tion between the spins, separated by a single qubit ro- 
tation. In present quantum dot systems, such a combi- 
nation of operations turns out to be rather difficult due 
to the strong influence of nuclear spins on the electron 
spin combined with the rather weak coupling between 
electron spins and electromagnetic fields. A more direct 
and robust C-Phase gate realization would therefore be 
desirable. 

New possibilities for constructing efficient two-qubit 
gates between single-spins are offered [§] by the recently 
demonstrated possibility to engineer different Zeeman 
splittings in neighbouring dots. Using micromagnets 
placed closed to the quantum dots, up to a 30 mT dif- 
ference in the local magnetic field in neighbouring dots 
has been demonstrated [9l-[llj|. and a 50 mT difference 
field is realistic 12]. By local control of the nuclear field, 



even higher difference fields, up to 1 Tesla, have been 
been achieved [ljj]. Also differences in g-factor between 
neighbouring dots offer similar possibilities and have been 
realized, both in quantum dots formed in etched GaAs 
pillars [14J and in InAs nanowires 15j. The main moti- 



vations for these gradients in Zeeman splitting were local 
addressing of individual spins by spectral selectivity, and 
the possibility of single-spin rotations by electric-dipole 
spin resonance [13]. Strong gradients also permit very 
fast exchange-controlled single-spin rotations 16]. 

Here we show that a difference in the local Zeeman 
splitting dE z between two quantum dots can be exploited 
to realize a fast and high-fidelity C-Phase gate under re- 
alistic conditions. Taking a broader view, we show that 
the natural two-qubit gate in tunnel-coupled quantum 
dots evolves from the exchange gate at dE z = to the 
C-Phase gate at large dE z . The cross-over depends on 
the strength of the interdot tunnel coupling and on the 
magnitude and rise time of the pulse that tilts the double 
dot potential. We quantify the error probability of the 
C-Phase gate under realistic conditions, giving insight 
in the optimal working point for fault-tolerant quantum 
computation. 

At the heart of the C-Phase gate are spin dependent 
energy shifts, analogous to C-Phase gate realizations with 
trapped ions [13] and nuclear spins in molecules [3]. As 
we will show in detail below, in the presence of a gradient, 
tilting the double dot potential lowers the two-electron 
states with anti-parallel spins in energy with respect to 
those with parallel spins. This results in energy shifts 
AE^i and AE±f respectively for the two antiparallel spin 
states. Fast and accurate control of the double dot poten- 
tial therefore allows adding a controllable spin-dependent 
phase shift to the anti-parallel spin states, described by 
the unitary transformation 
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expressed in the {| ft), I t-l-), I 4-t), I 4-4-)} basis and 
0i(2) = AE n ^t w(nt /h where t wait is the time the sys- 
tem spent in the tilted position. 

For cf>i + tp2 = 7T, the gate corresponds to a controlled 
7r-phase gate, up to single-qubit z rotations, with phase 
<f>\ and <j>2 on the first and second qubit respectively. The 
single-qubit z rotations can be left out when the C-Phase 
is part of a larger sequence of operations, as they can be 
absorbed into the phase of the subsequent single-spin x 
or y rotations (lif . 

In our analysis, we consider a closed double quantum 
dot filled with two electrons. The relevant two-electron 
states are the four states with one f or I electron in 
each dot and the singlet states 5(0, 2) (5(2, 0)) with both 
electrons in the right (left) dot. We also define 5(1, 1) = 
(I U) - I ;t»/V2 and T (l,l) = (| U) + I 4-t))/V2. 
The triplet states with two electrons in the same dot are 
much higher in energy and can be ignored in the following 
discussion. The Hamiltonian of this system in the basis 

{r+ = in), in), i n), r_ = \u), s(o,2), 5(2,0)} is 

then given by 
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where t is the interdot tunnel coupling (at the anti- 
crossing, the 5(1, l)-5(0, 2) separation is 2^/2t), E z is 
the average and dE z the difference in the Zeeman split- 
ting between the two dots, U the energy cost for moving 
both electrons into the same dot (analogous to the on- 
site interaction energy in the Hubbard model [Hf) and e 
the detuning or relative alignment of the potential of the 
two dots. In this Hamiltonian, we neglect the influence of 
spin-orbit interaction and magnetic field gradients in the 
x and y direction. The effect of the transverse gradients 
is to couple T + and T_ to 5(1,1), and indirectly (and 
weakly) to 5(0, 2) as well. The spin-orbit interaction 
gives a direct matrix element between the (1,1) triplet 
states and 5(0,2) [2(J. In both cases, the resulting T±- 
5(0, 2) anti-crossing will only produce a (small) energy 
shift and phase evolution of T±. Since this phase shift 
will be different from the phase shift of the anti-parallel 
spin states, we still obtain a C-Phase gate. 

The resulting energy level diagram as a function of 
detuning e is depicted in Fig. QJi. For e close to zero, 
the electron spins in the two dots are uncoupled and the 
eigenstates are T+, | t4-)> I 4-t) an d T— Their energy 
differences are set by E z and dE z . The T + and T_ states 
don't move in energy with e, since they are not coupled 
(via tunneling) to S(0,2) or S(2,0). The anti-parallel spin 
states, in contrast, shift down in energy as |e| moves away 
from zero, since they contain a 5(1, 1) component and 
the 5(1, 1) state has an avoided crossing with the 5(0, 2) 
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FIG. 1: (a) Energy level diagram for a double quantum dot 
containing two electrons, as a function of detuning e along the 
(2, 0)-(l, l)-(0, 2) axis. The six levels with the lowest energies 
are shown, and are labelled in the figure. Energies are ex- 
pressed relative to the (1, 1) states at e = and E z = 0. The 
parameters used are U — 1000/ieV, t — 10peV, dE z — 14/xeV 
and E z — 100/ieV. (b) Energy shift of the levels around e = 
with respect to their energy at e = 0. (Inset) Schematic 
double-dot charge stability diagram around the (1, 1) charge 
configuration. The dashed line shows the detuning axis along 
which panel (a) is plotted. 



and 5(2, 0) states. By controlling the detuning e in time, 
one can add a phase shift between the states with anti- 
parallel and parallel spins and therefore realize in a single 
step the C-phase gate presented in Eq. [TJ 

In the following, we will quantify the spin-dependent 
phase shifts from an effective Hamiltonian acting on the 
four two-qubit states. Analytical formulas for the phase 
accumulation <f>i and </>2 in Eq.Q]will help to characterize 
the speed and fidelity of the C-Phase gate. 

Using a Schrieffer-Wolf transformation and consider- 
ing t smaller than U ± e (> 0) [U |22| |. we can derive 
an effective Hamiltonian H e = e s He~ s eliminating at 
any order of approximation the off-diagonal terms of H 
between the (1,1) qubit states and 5(0,2) and 5(2,0). 
To second order, and expressed in the same basis as H 
in Eq. [2l 5 has the form 
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where j(dE z ) = t/(U - e- dE z /2) and a(dE z ) =t/{U^ 
e—dE z /2). Then, the projection of H e onto the two-qubit 
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FIG. 2: (a) Probability of | ft) i n one of the anti-parallel 
eigenstates which shows the change in the character of the 
eigenstates when dE z is increased, extracted from diagonali- 
sation of H e . Different curves correspond to values of e — 0, 
100, 750, 1200, 1600, 1800 (in red), and 1900 (U = 2000/xeV, 
t = 5fieV). The closer e is to U, the larger the dEz that 
is needed to change the eigenbasis. (b) Energy shift of the 
two anti-parallel qubit states as a function of dE z (extracted 
from Eq. [7} when e is pulsed from to e* = 1800/ieV 
(U — e* — 200/ieV, corresponding to the red line in the top 
panel) . For small dE z , the eigenstates are triplet and sin- 
glet and only one state (the singlet) is shifted in energy. At 
dE z — y2t, the energy shift is the same for the two anti- 
parallel eigenstates. 
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In the limit where e, dE z <C U, a and /3 are well approx- 
imated by 2t 2 /U and Eq. H reduces to Eq. 37 of Ref. 
with exchange interaction J = 4t 2 /U. 

When dE z is non-zero but small with respect to t, the 
energy eigenstates are { | tt) > I It) } close to e = but 
{S(l, 1), T (l, 1)} for larger e. This can be seen in Fig. Hi, 
where the tl component of one of the anti-parallel eigen- 
states is shown as a function of dE z for 7 values of e. The 



exchange gate takes advantage of this condition by puls- 
ing e non-adiabatically from zero to a finite value 0. 
Then constructing a C-Phase gate requires combining 
two exchange gates with a single-spin rotation. 

When dE z is large, the magnetic field gradient is the 
dominant part of the Hamiltonian, hence {| tl)J It)} are 
eigenstates as long as e is smaller than U, as seen also in 
Fig. [2jt. In this case, only the energies of the antiparallel 
states change. The exchange gate Qis then no longer 
possible and the corresponding eigenenergies vary with 
e, as seen from Fig. [5]d and from the following analytical 
expressions: 

a(dE z ) + a(-dE z ) , , s 

E nm (dE z ,e) = y ' 2 y T A(dE z ) (7) 

A(dE z ,e) = ll2jdEl + dE z Z(dE z ,e) + n{dE z ,e) (8) 



where T,(dE z ,e) = 2(a(—dE z ,e) — a(dE z ,e)) and 
Tl(dE z ,e) = 2{a(dE z ,e) 2 + a{-dE z ,e) 2 ). 

This means that a C-Phase gate can be realized in a 
single step via a fast pulse in e, from e = to e = e*. 
The closer to the 5(1, l)-5(0, 2) crossing the accumu- 
lation of phase takes place, the larger the energy shift 
hence the faster the operation is. We note that when 
dE z w \[2t, A — dE z to second order, so | ti) an d 
| |t) acquire the same energy difference A-E-f-m-M = 
E\\.(i-\){dE z , e*) — E^^Y[)(dE z ,§) and thus accumulate 
the same phase when e is pulsed. In this case, the C- 
Phase gate is obtained for <j>\ = 2 = tt/2 (as noted be- 
fore, the general condition for C-Phase is <j>\ + <j>2 = t). 
A ten ns gate duration, requires an energy shift of ~ 0.1 
/j,eV, which is in reach with dE z = 5/j,eV as seen from 
Fig. |3] (we note that phase evolution during the finite rise 
time of the pulses can be easily included in the timing 
calibration). Gate durations of 1 ns and less are pos- 
sible by using stronger tunnel couplings or larger pulse 
amplitudes. If dE z is made larger too, this will not af- 
fect the gate operation. These gate durations are 10 4 -10 5 
times shorter than reported spin-echo times in GaAs [f| , 
so error probabilities due to decoherence will be below 
the 10 -4 accuracy threshold for fault-tolerant quantum 
computation 

Finally, we analyze for realistic conditions the gate er- 
rors that can be expected from other sources than deco- 
herence, with the help of the analytical expressions ob- 
tained above. To simplify the analysis, we consider the 
case where <f>\ = 4>2 = <\>. As in the case of the exchange 
gate 4> will be sensitive to potential fluctuations due 
to charge noise or gate voltage noise since we allow a 
small amount of mixing between charge and spin states 
to obtain the phase shift. The effect of potential fluctua- 
tions manifests itself predominantly as fluctuations in e, 
rather than as fluctuations in t or U 25]. The closer to 
the 5(1, l)-5(0, 2) crossing the operating point is located, 
the more sensitive to detuning fluctuations the C-Phase 
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FIG. 3: Expected energy shift with respect to the case e = 
as a function of detuning e* (t = dEz/V% = 5/ieV and 
U = 2000^teV). The closer to the 5(1, l)-5(0, 2) anti-crossing, 
the faster the gate. Also shown is the relative errors in the 
phase <f> due to lfieV fluctuations in e. Between e* = 300/ieV 
and e* = 1850/ieV, errors are below the threshold for fault- 
tolerant quantum computation. 



gate is, as can be seen from Fig. El From Eq.0 analytical 
expressions can be derived for d/fi/cf) induced by detuning 
fluctuations and the results are also shown in Fig. |3l tak- 
ing a realistic value of IfieY detuning fluctuations [26j j . 
We see that relative errors dtfi/cj) down to 6 x 10~ 3 can 
be obtained for the 10 ns C-Phase gate, which translate 
to a gate error probability (d<p/4>) 2 = 4x 1CP 5 . 

A second possible source of errors is the random value 
of the (quasi-static) nuclear field in the two dots. This 
affects the time-evolution during C-Phase in two ways. 
First, it causes a dephasing that is independent of e. This 
dephasing has been shown to be (largely) reversible by 
spin-echo techniques and since C-Phase commutes with 
dephasing, it can be simply embedded in existing spin- 
echo sequences. Second, there is an additional contribu- 
tion to the e-dependent energy shift <f>, but this is only 
a second order effect in the dE z fluctuations. The slow 
variation of <f) with dE z is seen in Fig. [5J in particular as 
dE z approaches t. Gate error probabilities due to dE z 
fluctuations can be evaluated from Eq. [7] in the regime 
of fast gate operation where dE z <C U — e* <C U and, for 
dE z = y/2t, are approximately given by d<j>/<j) = 6/ (U—e) 
where S is the size of the fluctuations in dE z (ss 20 neV 
in GaAs). Therefore, for the values used above to obtain 
a 10 ns gate, the gate error probability (d<p/4>) 2 due to 
the nuclear field randomness would be only ss 10~ 8 . 

The last source of error we consider is possible non- 
adiabaticities that may occur if the rise time of the de- 
tuning pulses is too short. This is relevant when the anti- 
parallel eigenstates are slightly different between the two 
pulse positions e = and e = e*. For the parameters 
used above, we quantify this difference via the square of 
the overlap between the eigenstates at the two pulse po- 
sitions and it is different from 1 only by 8 x 10~ 4 . For 
instantaneous detuning pulses, the resulting gate error 
probability calculated from simulations is of the same 
order of magnitude, ~ 2 x 10~ 3 . However, for realistic 



rise times in the range of 0.3-1 ns and again using the 
same t and e as above, the gate error probability is well 
below 10~ 4 . By using shaped pulses, the adiabaticity can 
be improved even further (27^ . 

In conclusion, we present and analyze a one-step im- 
plementation of a C-Phase gate between two single-spin 
qubits in a double quantum dot. This gate is the nat- 
ural replacement of the exchange gate when large mag- 
netic field or ^-factor gradients are present. We estimate 
that gate operation times of a few nanoseconds can be 
achieved under realistic conditions with fidelities above 
the accuracy threshold for fault-tolerance. To probe the 
spin-dependent phase shifts produced by the C-Phase 
gate, a Ramsey type experiment can be used, with the 
C-Phase gate in between two single-spin rotations. We 
expect that this work will impact ongoing experiments 
aimed at the realization of elementary quantum gates 
and simple quantum protocols with single-spin qubits in 
a variety of quantum dot systems. 
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